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ABSTRACT 

We present deep Chandra observations and Spitzer Space Telescope infrared (IR) spectroscopy of 
the shell in the composite supernova remnant (SNR) Kes 75 (G29.7-0.3). The remnant is composed 
of a central pulsar wind nebula and a bright partial shell in the south that is visible at radio, IR, 
and X-ray wavelengths. The X-ray emission can be modeled by either a single thermal component 
with a temperature of ~ 1.5 keV, or with two thermal components with temperatures of 1.5 and 0.2 
keV. Previous studies suggest that the hot component may originate from reverse-shocked SN ejecta. 
However, our new analysis shows no definitive evidence for enhanced abundances of Si, S, Ar, Mg, 
and Fe, as expected from supernova (SN) ejecta, or for the IR spectral signatures characteristic of 
confirmed SN condensed dust, thus favoring a circumstellar or interstellar origin for the X-ray and IR 
emission. The X-ray and IR emission in the shell are spatially correlated, suggesting that the dust 
particles are collisionally heated by the X-ray emitting gas. The IR spectrum of the shell is dominated 
by continuum emission from dust with little, or no line emission. Modeling the IR spectrum shows 
that the dust is heated to a temperature of ~ 140 K by a relatively dense, hot plasma, that also gives 
rise to the hot X-ray emission component. The density inferred from the IR emission is significantly 
higher than the density inferred from the X-ray models, suggesting a low filling factor for this X-ray 
emitting gas. The total mass of the warm dust, component is at least 1.3 x 10“ 2 Mg, assuming no 
significant dust destruction has occurred in the shell. The IR data also reveal the presence of an 
additional plasma component with a cooler temperature, consistent with the 0.2 keV gas component. 
Our IR analysis therefore provides an independent verification of the cooler component of the X-ray 
emission. The complementary analyses of the X-ray and IR emission provide quantitative estimates 
of density and filling factors of the clumpy medium swept up by the SNR. 

Subject headings: dust, extinction - infrared: ISM - ISM: individual objects (SNR G29.7-0.3) - ISM: 
supernova remnants - pulsars: individual (PSR J1846-0258) - X-rays: ISM 


1. INTRODUCTION 

The efficiency of dust formation in supernovae (SNe) 
and its subsequent evolution are still not well under- 
stood. Large discrepancies exist between theoretical 
predictions for dust formation in SNe and dust masses 
measured from recent infrared (IR) observations of su- 
pernova remnants (SNRs). Theoretical models based 
on the classical nucleation theory (Nozawa et al. 2003; 
Kozasa et al. 2009) and the chemical kinematic approach 
for the formation of molecular precursors to dust in SN 
ejecta (Cherchneff & Dwek 2010) predict that the total 
dust mass produced per SN explosion should be 0. 1—0.7 
M©. A large fraction of this dust may be destroyed 
by erosion and sputtering between the forward and re- 
verse shocks, but the surviving mass is still expected to 
be 0.04-0.2 M© (e.g., Dwek et al. 2008; Kozasa et al. 

2009) . In addition, dust evolution models show that 
an average SN would need to produce 0. 1-1.0 M© of 
dust in order to explain the large quantities of dust ob- 
served in high-redshift galaxies (e.g., Michalowski et al. 

2010) . However, even such large production efficien- 
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cies may not be sufficient to explain the mass of dust 
in these objects (Dwek & Cherchneff 2011; Dwek et al. 
2011; Valiante et al. 2011; Gall et al. 2011a, b). 

In recent years, observations of SNRs with Spitzer 
have provided mass estimates of newly formed warm 
(70-200 K) dust for several remnants, but these esti- 
mates lie in the 0.02-0.1 M© range (e.g., Sugerman et al. 
2006; Rho et al. 2008, 2009; Temim et al. 2010), signif- 
icantly less than theoretical predictions. Recent Her- 
schel observations revealed 0.4-0. 7 M© of cool dust in 
SN 1987A, the highest mass estimate for freshly formed 
dust up to date (Matsuura et al. 2011). Spitzer, Her- 
schel and AKARI observations of Cas A imply that 
0.08 M© of dust is present in the remnant (Rho et al. 
2008; Sibthorpe et al. 2010; Barlow et al. 2010), consis- 
tent with model calculations for the evolution of dust 
formed in the ejecta of Type lib SNe (Nozawa et al. 
2010). Spitzer observations of G54. 1+0.3 showed evi- 
dence for up to 0.1 M© of fresh dust surrounding the 
pulsar wind nebula (PWN; Temim et al. 2010). How- 
ever, a large fraction of this dust may be destroyed by 
the SNR, reverse shock, so the total mass estimate still 
falls short of what theories predict. 

In order to determine the net amount of dust injected 
by SNe into the interstellar medium (ISM), it is also im- 
portant to understand the dust destruction rates in SNR 
shocks. Hydrodynamic simulations of dust destruction 
in SNR reverse shocks suggest that 20%-100% of dust is 
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Fig. 1. Three-color image of Kes 75, where the Chandra X-ray 
emission is shown in blue, MIPS 24 /cm emission in green, and radio 
emission in red. The IR emission shows a partial shell in the SE 
and SW regions and coincides with the X-ray and radio emission. 
The saturated point sources seen in green are foreground stars. 


destroyed by the shock, with the fraction depending on 
the grain species, ambient gas density, and the type of SN 
(Nozawa, et al. 2007; Silvia et al. 2010). The simulations 
show that grains smaller than 0.1 pm are either sputtered 
to smaller radii or completely destroyed. Recent model- 
ing of the IR emission from Puppis A and the Cygnus 
Loop shows clear evidence of dust grain destruction in 
an SNR shock (Arendt et al. 2010; Sankrit et al. 2010). 
The study of Large Magellanic Cloud (LMC) SNRs for 
both Type la and core-collapse SNRs by Borkowski et al. 
(2006) and Williams et al. (2006) shows that 40% of the 
total mass in dust grains has been destroyed, including 
as much as 90% of all grains smaller than 0.04 pm. The 
inferred dust-to-gas mass ratios for these SNRs were sev- 
eral times smaller than the typically assumed LMC value 
of 0.025. In order to explain these results, either a much 
higher dust destruction rate or lower preshock dust-to- 
gas mass ratios is required (Borkowski et al. 2006). 

Kes 75 is an important subject for multi-wavelength 
studies of the dust and gas content of SNRs. Because of 
its relatively young age, it allows the possible distinction 
between ejecta and swept-up material giving rise to the 
X-ray and IR emission. Furthermore, it enables studies 
of the shocked dust and gas properties and grain destruc- 
tion efficiencies in shocks. In a recent study of Chandra 
and Spitzer observations of Kes 75, Morton et al. (2007) 
suggest that the X-ray emission in Kes 75 arises from 
reverse-shocked ejecta, and also estimated the dust-to- 
gas mass ratio for the SNR shell to be ~ 10 -3 , several 
times lower than what is typical for the Galaxy. However, 
these estimates were based only on the Multiband Imag- 
ing Photometer (MIPS) 24 pm measurement. Kes 75 was 
not detected at Infrared Array Camera (IRAC) wave- 
lengths (3-8 pm), suggesting that the small grains in the 
shell have been destroyed by shocks (Morton et al. 2007). 
Our new Spitzer spectroscopy allows us to place more 
detailed constraints on the dust properties and mass by 
fitting the spectra from the shell directly. 


Kes 75 is a Galactic SNR consisting of a PWN and a 
partial thermal shell. The PWN is powered by a 0.3 s 
pulsar (PSR .11846-0258) that has a characteristic age 
of 723 yr, one of the youngest known (Gotthelf et al. 
2000). An upper limit on its true age of 884 yr was 
derived from the measurement of the braking index, 
n = 2.65 ±0.01 (Livingstone et al. 2006). The spin-down 
rate P — 7.1 x 10~ 12 ss -1 leads to a spin-down luminosity 
of E = 8.1 x 10 36 erg s _1 and a strong surface magnetic 
field B = 5 x 10 13 G (Gotthelf et al. 2000). The dis- 
tance to Kes 75 is poorly constrained, with the estimates 
ranging from 5.1 to 21 kpc (Caswell et al. 1975; Milne 
1979: Becker & Helfand 1984; Leahy & Tian 2008). Re- 
cent millimeter observations of CO line emission toward 
Kes 75 provide evidence for an association with an adja- 
cent molecular cloud at a distance of ~ 10.6kpc (Su et al. 
2009). We adopt this distance throughout the paper, at 
which the radius of the SNR shell is equal to 5.7 pc. The 
dependence of the derived gas and dust properties on the 
assumed distance will be discussed later in the paper. 

Radio observations of Kes 75 show a partial SNR shell, 
~ 1(5 in radius, with a spectral index a = 0.7, and a flat- 
ter PWN component with a = 0.25 (Becker & Kundu 
1976; Helfand et al. 2003). Here, a is defined such that 
L v ex v~ a , where L„ is the synchrotron luminosity. X- 
ray observations with A SC A and Chandra also showed 
the composite nature of the remnant with a morphology 
closely resembling the radio (Blanton & Helfand 1996; 
Helfand et al. 2003). Chandra spectroscopy of the PWN 
shows that the emission softens with distance from the 
pulsar, with the photon index, F ranging from 1.1 to 
1.9 (Morton et al. 2007; Ng et al. 2008), where F re- 
lates the number of photons N and energy E such that 
dN/dE c x E~ r . The shell emission in Kes 75 is primarily 
concentrated in two regions in the southeast (SE) and 
the southwest (SW). Chandra observations show that 
the spectra of the shell can be characterized by a two- 
temperature thermal model, possibly associated with the 
forward-shocked material and the reverse-shocked ejecta 
(Morton et al. 2007). 

The first IR detection of Kes 75 with Spitzer revealed 
a partial IR shell coincident with the shell at radio and 
X-ray wavelengths (Morton et al. 2007). Morton et al. 
(2007) used Chandra and MIPS 24 pm observations, 
along with IRAC upper limits, to constrain the dust 
emission models and obtain a dust-to-gas mass ratio 
Mdust/Mgas ~ 5 x 10~ 4 , more than an order of mag- 
nitude lower than that of the Galaxy. Since Kes 75 was 
only detected in the MIPS 24 pm band, the data were not 
sufficient to characterize the dust composition and place 
tight constraints on the dust temperature and mass. 

Here, we present Chandra X-ray observations and 
Spitzer IR spectroscopy of the shell in Kes 75 that al- 
lowed us to better characterize the observed dust and 
gas properties. Fitting of the X-ray spectra allows us to 
determine the properties of the thermal emission from 
the shell, while the modeling of the IR spectra from the 
shell provides the dust mass and temperature, and an 
independent estimate of the temperature and density of 
the heating plasma. The comparison of the X-ray and 
IR modeling results provides additional insight into the 
nature of the shell emission in Kes 75. 

2. OBSERVATIONS AND DATA REDUCTION 
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TABLE 1 

MIPS 24 MICRON Fluxes 


Region 

Flux (Jy) 

Ext. Corr. Flux (Jy) 

SE Rim 

1.9 ± 0.4 

5.1 ± 1.0 

SW Rim 

2.7 ± 0.5 

7.0 ± 1.4 

Emission from the Shell minus Inner SNR 

SE Rim 

0.9 ± 0.2 

2.2 ± 0.4 

SW Rim 

1.3 ± 0.3 

3.4 ± 0.7 


Note. — The uncertainties include MIPS 
calibration uncertainties, but do not account, 
for the uncertainties in the extinction correc- 
tion. The extinction correction was applied 
using the extinction curve of Chiar & Tielens 
(2006) and a hydrogen column density of Nfj = 

3.5 x 10 22 cm -2 . 

2.1. Chandra 

The Chandra observations of Kes 75 were carried out 
with the Advanced CCD Imaging Spectrometer on 2006 
June 5, 7, 9, and 12, under observation II) numbers 7337, 
6686, 7338, and 7339, and an exposure time of 18, 55, 40, 
and 42 ks, respectively. The dataset was previously an- 
alyzed by Morton et al. (2007) and Su et al. (2009), but 
here, we reanalyze the spectra using a different approach 
to background selection and fitting. The data were re- 
processed and cleaned using CIAO v. 4.3, resulting in 
total clean exposure time of 152 ks. 

Since the observations were carried out close in time 
and have similar chip orientation, we merged the indi- 
vidual data sets into a single event file. The merged 
image is shown in Figures 1 and 2. We chose extrac- 
tion regions for the SE and SW rims similar to those in 
Morton et al. (2007) and Su et al. (2009), with an addi- 
tional region covering the inner SNR between the PWN 
and the southeastern rim of the shell. The extraction 
regions for the X-ray spectra are shown in Figure 2(a). 

2.2. Spitzer Spectroscopy 

Spitzer IR spectroscopy of the shell in Kes 75 was 
carried out on 2008, June 5 (program II) 50447) using 
the long-low (LL) and long-high (LH) modules of the In- 
frared Spectrograph (IRS: Houck et al. 2004) that cover 
the range from 14-37 pm. The locations of the IRS 
slits overlaid on the MIPS 24 pm image (Morton et al. 
2007) are shown in Figure 3. The spectra were taken 
at positions 18 ! )46i n 22!3, -2°59'46"9 and lStMei^OiSe, - 
2 0 58'54"9, corresponding to the SW and SE portions of 
the shell, respectively. In order to measure the back- 
ground emission for the LH module, we also placed a slit 
at position 18 1 )46( n 24*03. -3°01 l '49'/0. south of the partial 
shell. All observations were taken with 10 cycles of 6 s 
exposures and the data were processed with the pipeline 
version SI 8. 0.2. 

The LH spectroscopy was cleaned with the IRS 
Rogue Pixel Mask Editing and Image Cleaning Soft- 
ware (IRSCLEAN1.9). The background spectrum was 
subtracted from the LII data and the Basic Calibrated 
Data of each nod position were combined. Spectra for 
the SW and SE portion of the shell were extracted us- 
ing the Spitzer IRS Custom Extractor (SPICE) with full 
slit extractions and extended source calibration. For the 



Fig. 2. (a) Three-color Chandra X-ray image of Kes 75, with 

0.3-2.0 keV emission in red 2.0-4.0 keV emission in green and 
4. 0-9.0 keV emission in blue. The extraction regions for X-ray 
source and background spectra are shown in white, (b) Two-color 
images of Kes 75 with MIPS 24 pm emission in red and Chandra 
X-ray emission from 0.5— 9.0 keV in cyan. The sub-slits of the IRS 
LL module used for the low-resolution IR spectral extraction are 
shown in white. We note that the regions labeled “Inner” in top 
and bottom panels do not spatially coincide. 

LL modules, the spectra were extracted using the CUbe 
Builder for IRS Spectral Maps (CUBISM; Smith et al. 
2007). which allowed us to select sub-slits from the off- 
source nods for background subtraction. Since the back- 
ground in the vicinity of Kes 75 is spatially variable 
(89-102 MJy sr -1 at 24 pm), some off-source nod po- 
sitions produced negative residuals in the background- 
subtracted spectra. The final background was produced 
by averaging spectra extracted from off-source sub-slits 
that did not result in negative residuals. The source 
spectra were extracted from sub-slits covering the shell 
emission and inner region on the SNR. and are shown 
in Figure 2(b). The rest of the analysis was carried out 
with the Spectroscopic Modeling, Analysis, and Reduc- 
tion Tool (SMART; Higdon et al. 2004), including fitting 
of the emission lines. 
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Fig. 3. Positions of the IRS slits overlaid on the MIPS 24 pm 
image. The LL slits are the large slits oriented in the east/west. 
direction and span the entire length of the SNR shell. The LH 
slits are the smaller slits shown in yellow, positioned at the east 
and west portions of the shell. The LH slits located south of the 
remnant were used for background subtraction. 

3. ANALYSIS AND RESULTS 
3.1. Multi-wavelength Morphology 

The three-color image of Kes 75, displayed in Figure 
1, shows the Chandra X-ray emission in blue, MIPS 24 
p m IR emission in green, and Very Large Array radio 
emission in red (Becker & Helfand 1984). The IR image 
shows a partial SNR shell, with the brightest emission 
concentrated in the SE and SW rims which coincides 
spatially with the X-ray and radio emission. The 24 
pm image also shows fainter IR emission in the inner 
region of the SNR, between the PWN and the south- 
ern rim, that is more diffuse than the filamentary shell. 
Extinction-corrected 24 pm fluxes of the bright portions 
of the shell are measured to be 5.1 and 7.0 Jy in the SE 
and SW rims, respectively. We also measured the total 
flux in the shell using the diffuse IR emission inside of the 
SNR as background. Since the IR spectra suggest that 
the emission from the shell is composed of two distinct 
dust temperature components, the latter flux densities 
are estimates for the warm dust component that, is con- 
centrated in the shell (see Section 3.3). The measured 
values are listed in Table 1. As discussed in Section 3.3, 
the IR emission is dominated by dust continuum, with 
very little or no contribution from line emission. The 
radio synchrotron emission from the shell has a spectral 
index a — 0.7 and a flux density of ~ 10 Jy at 1 GHz 
(Becker & Kundu 1976). When extrapolated to Spitzer 
IR wavelengths, the contribution from the synchrotron 
emission is negligible. 

The X-ray emission spatially correlates with the IR, as 
is best seen in Figure 2. The rims of the SNR shell are 
dominated by thermal emission with a range of temper- 
atures. This is best seen in the X-ray three-color image 
in Figure 2(a). The colors show that some portions of 
the shell are dominated by softer X-ray emission, while 


other regions have a harder X-ray spectrum (blue). The 
SE portion of the shell appears particularly clumpy, with 
softer X-ray clumps embedded in more diffuse hard emis- 
sion. Figure 2 also shows that diffuse X-ray emission fills 
the inner SNR. Su et al. (2009) show that the spectrum 
of the diffuse emission north of the PWN is featureless 
and most likely non-thermal in nature. This emission 
likely represents the dust-scattered halo from the PWN. 
However, the emission from the inner SNR region south 
of the PWN is more enhanced than in the north, and 
shows evidence of thermal X-ray emission. 

3.2. X-Ray Spectroscopy 

The Chandra X-ray spectra were extracted from re- 
gions shown in Figure 2(a). The best fit parameters 
are shown in Table 2. We also refitted the spectrum of 
the entire PWN, excluding the pulsar, with a power-law 
model and the Tuebingen-Boulder ISM absorption model 
XSTBABS. The spectrum contained ~ 100,000 counts 
after subtracting a background outside of the SNR,. The 
absorbed power-law model provided a good fit to the 
spectrum, giving a y 2 of 0.92. Using the XSTBABS ab- 
sorption model, we find Nr — 3.75 x 10 22 em -2 , a slightly 
lower value than reported previously, and F = 1.98±0.3, 
consistent with previous results. 

Su et al. (2009) modeled the X-ray spectra from the 
SE and SW rims with a thermal component, plus a non- 
thermal one to account for possible dust scattering. In- 
stead of using this approach, we assumed that the emis- 
sion from the dust scattering halo is approximately sym- 
metric about the PWN, and we chose background regions 
in the northern part of the shell that do not appear to be 
contaminated by thermal emission to subtract the scat- 
tered emission component. The background regions are 
shown in Figure 2(a). The region labeled “Bkgl” was 
used as background for the “Inner” region, while the re- 
gion labeled “Bkg2” was used as background for the SE 
and SW regions. The source spectra for the SE (SW) rim 
contained 37,000 (16,000) counts, while the total counts 
in the background, including emission from the dust scat- 
tering halo, was approximately 2,900 (2700) counts. The 
total number of source counts in the “Inner” region was 
8900, but the background in this case was significantly 
higher, approximately 5350 counts. The spectra from the 
SE and SW regions were grouped by 30 counts bin -1 , 
while the spectrum from the inner region was grouped 
by 100 counts bin -1 . 

The X-ray spectra were fit with two different mod- 
els, a non-equilibrium ionization collisional plasma model 
(VNEI) and a plane-parallel shock plasma model (VP- 
SHOCK), that also include elemental abundances as free 
parameters. The models are shown in Figures 4, 5, 6, 
and 7. A single temperature VNEI thermal component 
did not provide a good fit to the SE and SW regions, 
while these areas were well-fit by the single-temperature 
VPSHOCK model (see Figures 4 and 7). The best-fit pa- 
rameters for all models are listed in Table 2. The single- 
temperature VPSHOCK model implies a temperature of 
1.48 keV and 1.93 keV in the SE and SW, respectively, 
and 2.47 keV in the inner SNR region. The ionization 
timescales are on the order of 10 n scm -3 in the shell and 
i0 10 sem -3 in the inner SNR region. We let the hydrogen 
column density vary and found that it is somewhat lower 
than for the PWN {Nr — 3 x 10 22 cm -2 ), as suggested 
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TABLE 2 

Chandra Spectral Fitting Results 


Parameter 

PUN 

SE Shell 

SVV Shell 

Inner SNR 


VPSHOCK-One-component Model 

jV//(10 22 cm“ 2 ) 

q 'tc+ 0-06 
o, '°-0.06 

o qo+0.02 

z -y«J_0.04 

q I q+ 0.03 
°- iZ -0 .10 

q tro+0.22 

*>'Oo_ 0 _23 

Photon Index (P) 
kT(keV) 

r.yo_ 0 .03 

1 4&40-02 

i -4o_ 0 02 

i-93tS;g 

o 47+ 0 - 22 

Z - 4 '~0.23 

t(s cm -3 ) 


(1.46!°;^) x 10 n 

(7.36! 3 ; 2 ^) x 10 10 

(1.65±g;||) x 10 10 

F(unabs. flux (erg cm -2 s -1 )) 

3.6 x 10” u 

5.5 x 10~ u 

3.2 x 10” 11 

1.8 x 10”” 11 

Normalization 


(7.7«; 3 ) x 10”” 3 

(2.8i±g;§|) x 10- 3 

(8-1 til) * 

Reduced x 2 statistic 

0.92 

1.33 

1.34 

0.92 

VPSHOCK-Two-component Model 

A(f/(10 22 cm 2 ) 


q oq+0. 18 

4 -°oig;i7 


fcTj (keV) 


i-46tg;S 



ri (s cm -3 ) 


(1.09!“;!®) x 1QU 

(5.23^;“) xlO 10 


Fi(unabs. flux (erg cm -2 s -1 )) 


7.23 X 10- 11 

6.1 x 10- 11 


Normalization! 


(8.6t°; 3 ) x 10 " 3 

(4.6ig;J) X 10-3 


fcT 2 (keV) 


°-20tg;S 

0 14+ 002 
u,14 -0.01 


r 2 (s cm -3 ) 


(2-97 x 10 9 

(2.22lJ;J|) x 10 10 


F 2 (unabs. flux (erg cm -2 s -1 )) 


3.4 x 10” 9 

1.0 x 10“® 


Normalization 


2 - 8 tl‘.l 

14.3±|;| 


Reduced x 2 statistic 


1.01 

0.92 


VNEI Model 

ATjf(10 22 cm 2 ) 


q oo'i'0.02 
0.01 

3.49tr 0 ® 

q 7 q-f0.25 

,3 * ' * 3 — 0.24 

kT i(keV) 


i-eorggf 

2-08l®-j; 

Oil 'f0.55 
z,ii -0.37 

ri(s cm -3 ) 


(5-7otS:S) x 1010 

(3.14±g;||) x 10 10 

(7-06ti;|?) X 10® 

Fi(unabs. flux(erg cm" 2 s -1 )) 


4.0 x 10- 11 

2.8 x 10“ u 

2.6 x 10 ” u 

Normalization! 


(6.80±g;gf ) x 10- 3 

(2.5±g |) x 10“ 3 

(l.org j) x 10“3 

fcT 2 (keV) 


0 94+0-04 
U '^ 4 -0.01 

0 90+0-09 

u - zu ~-0 oQ 


T '2 (s cm - ' 3 ) 


> 1.5 x 10 13 

> 2.57 x 10 12 


F 2 (unabs. flux (erg cm” 3 s -1 )) 


2.1 x 10“ 10 

1.5 x ur 10 


Normalization 2 


0 141 +0-007 
U,i4i -0.005 

0 11+ 0 ' 6 


Reduced x 2 statistic 


1.41 

1.27 

0.97 


Note. — The listed uncertainties are 1.6 a (90% confidence) statistical uncertainties only. The apertures 
used in the spectral extractions are shown in Figure 2 (top), and the spectra and fits are shown in Figures 
4~6. The fluxes were calculated in the 0. 5-9.0 keV. The normalization of the thermal models is equal to 
1 0 ~ 1 4 n e ?iH V/47rd 2 cm" 5 , where V is the volume of the emitting region and d is the distance to the SNR. 
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Fig. 4. Chandra X-ray spectra and best-fit models for the SE (left) and SW (right) regions of the shell. The extraction regions are 
shown in Figure 2 and best-fit models listed in Table 2. The spectra were fitted with a single-temperature VPSHOCK model; 1.5 keV in 
the SE and 1.9 keV in the SW. Observed emission lines are from Si. S, Mg, Ar, and Fe. 



j : , . , . j 

1 10 

Energy (keV) 



j , . ■ ■ ■ ■ < 
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Energy (keV) 


Fig. 5. Same as Figure 4, but with a two-temperature VPSHOCK model (solid line), composed of a hotter component with a temperature 
of 1.5 (1.6) keV (dotted line) and a cooler component with a temperature of 0.2 (0.14) keV (dashed line), for the SE (SW) shell. 



Energy (keV) 



Energy (keV) 


Fig. 0. Same as Figure 4, but with a two-temperature VNEI model (solid line), composed of a hotter component with a temperature 
of 1.6 (2.1) keV (doited line) and a cooler component with a temperature of ~ 0.2 (0.2) keV (dashed line) for the SE (SW) shell. 
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Fig. 7. Chandra X-ray spectra and best-fit VPSHOCK model 
for the inner SNR region shown in Figure 2. The best-fit thermal 
temperature is ~2,5 keV and the observed emission lines are from 
Si, S, Mg, Ar, and Fe. Other parameters are listed in Table 2. 

by Su et al. (2009). 

As will be discussed in the following sections, the IR 
data imply the presence of an additional, lower tem- 
perature plasma component. A somewhat better \ 2 
is achieved using a two-temperature VPSHOCK model, 
shown in Figure 5. For comparison, we also fit the data 
with a two- temperature VNEI model (Figure 6). but the 
fit was slightly worse in this case. The temperature of 
this second component is approximately 0.2 keV for both 
models, but the normalization and ionization timescale 
are not well constrained. A similar y 2 value is achieved 
for a relatively short ionization timescale of ~ 10 9 scm -3 
and a very long timescale of ~ 10 13 s cm -3 . The F-test 
implies that the presence of this lower temperature com- 
ponent is not statistically significant. However, as will 
be discussed in the following sections, the IR data imply 
the presence of an additional, lower temperature plasma 
component. 

The X-ray spectra show clear evidence for Si, S, Ar, 
Mg, and Fe lines, but in contrast to previous studies 
(Morton et al. 2007; Su et al. 2009), we find no definitive 
evidence for enhanced abundances in either of the ther- 
mal components. The main difference in our analysis is 
a more careful approach in background selection that ac- 
counted for the emission from the dust scattering halo. 
After subtracting this emission from the source spectra, 
the enhanced Fe abundance was no longer required by 
the fit. While the best fit residuals in Figures 4, 5, and 
6 clearly show that the models do not provide a good 
fit to some of the emission lines, varying the abundances 
did not significantly improve the fit, and the abundances 
remained around the solar value. We believe that the 
residuals are caused by the variation in temperature in 
the higher energy component with a shorter ionization 
timescale. These differences in temperature are clearly 
seen in Figure 2(a), and indeed, fits to small sub-regions 
in the SE and SW shell show that the temperature ranges 
from 0.9 keV for some of the smaller clumps, up to 2 keV 
for more diffuse regions. The abundances were consistent 
with solar even for these various sub-regions. 

3.3. Infrared Spectroscopy 


The IRS slits were positioned along the brightest re- 
gions of the SW and SE parts of the shell (see Figure 
3). The final background subtracted LH spectra of the 
shell are shown in Figure 8. The spectra of the shell are 
dominated by continuum emission from dust with very 
little contribution from line emission. Emission lines that 
are present in the background-subtracted LII spectra are 
listed in Table 3, and they include [O IV] (25.89 pm), 
[Fe II] (25.99 pm), [S III] (33.48 pm), and [Si II] (34.82 
pm). The listed uncertainties in the observed line inten- 
sities are statistical uncertainties only, and they do not 
include the uncertainties caused by the spatially vari- 
able background in the vicinity of Kes 75. Since the 
background for the high-resolution spectra was extracted 
from a single region south of the SNR (see Figure 3), it 
is difficult to say how much line emission, if any, is as- 
sociated with Kes 75. For this reason, the measured line 
intensities should be treated as upper limits. All of the 
lines except O IV appear in the local background. In 
addition, the lines all have a full width at half-maximum 
(FWHM) that is similar to the expected resolution of the 
IRS of 500 kms' 1 , and none show evidence for signifi- 
cant broadening that would be expected from a rapidly 
expanding blast wave. There are no lines detected in 
the LL modules after background subtraction. However, 
the spectra have high uncertainties due to the high local 
background that varies significantly across the LL slits. 

Extinction correction was applied to the emission line 
intensities and the LL spectra used in the dust model 
fitting. We used the extinction curve of Chiar & Tielens 
(2006), and the relation Nh/Ak = 1.821 x 10 22 cm -2 
from Draine (1989). We set the value of A r n = 3.5 x 
10 22 cm" 2 , an average derived from the fitting of the 
Chandra data to the SE and SW shell. 

3.4. Modeling of the Dust Emission 

The morphological similarities between the X-ray and 
IR emission in the southern part of the shell suggest 
that the emitting dust is likely collisionally heated by 
the shocked, X-ray emitting gas. The lack of strong IR 
line emission also indicates that no cooler and denser 
regions are present where strong UV and optical emis- 
sion might radiatively heat the dust. The dust grain 
model used to fit the IR spectra is the BARE-GR-FG 
model from Zubko et al. (2004) that includes polycyclic 
aromatic hydrocarbons (PAHs), silicate, and graphite 
grains. Zubko et al. (2004) determined the grain size dis- 
tribution for the model by simultaneously fitting the av- 
erage interstellar extinction, emission, and abundances. 
Our fitting was performed using a dust heating code that 
calculates spectral models for dust grains immersed in 
an X-ray plasma (Dwek 1986: Arendt et al. 2010). The 
model includes stochastic heating of smaller grains and 
thermal sputtering of dust grains that reduces the size of 
the initial grain population. The free parameters of this 
dust model are the electron density n e , electron temper- 
ature T e , and the change in the size of each grain due to 
sputtering, Aa, that depends on the age of the shocked 
gas. We produced a grid of models using a density range 
n e (cm -3 )— [0.35, 0.7, 1.4, 2.8, 5.7, 11. 23, 45, 90, 180], 
and T e (K)=[1.0,1.8, 3.2, 5.6, 10, 18, 32jxl0 6 . The sput- 
tering amount Aa was varied from 0.0-0.01 pm in steps 
of 0.0005 pan. 
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TABLE 3 
IRS Line Fits 



Line ID 

Line Center 
(pm) 

Line Flux 

(10 ~ 6 erg cm 2 ~ s~ x sr~ 1 ) 

De-reddened Flux 
(I0~ 6 erg cm 2 ~ s~ l sr~ 1 ) 

FWHM 

(pm) 

FWHM 

(km s ~~ 1 ) 

LH east 







'O IV] (25.8903) 

25.8914 ± 0.0012 

13.47 ± 1.03 

31.52 ± 2.41 

0.047 ± 0.001 

544 


Fe II] (25.9883) 

25.9751 ± 0.0013 

2.86 ± 0.25 

6.66 ± 0.58 

0.044 ± 0.001 

511 


S Hit (33.4810) 

33.4733 ± 0.0019 

37.81 ± 4.51 

71.84 ± 8.57 

0.051 ± 0.002 

458 


Si II] (34.8152) 

34.8113 ± 0.0011 

37.13 ± 2.11 

68.80 ± 3.97 

0.056 ± 0.001 

482 

LH west 







[O IV] (25.8903) 

25.9009 ± 0.0019 

4.87 ± 0.58 

11.39 ± 1.36 

0.049 ± 0.002 

564 

[S III] (33.4810) 

33.4816 ± 0.0014 

24.06 ± 1.94 

45.71 ± 3.69 

0.054 ± 0.001 

483 


Note. — The listed uncertainties are ler statistical uncertainties from the fit only. Due the highly variable background 
in the vicinity of Kes 75, the measured line intensities should be treated as upper limits. 


We used the above models to fit spectra from three 
different regions along the Kes 75 shell; the brightest SE 
and SW rims and the fainter diffuse emission inside of the 
shell. The LL sub-slits used for spectral extractions are 
shown in Figure 2(b). We initially used a background 
from outside of the SNR for all three regions, and we 
found that the spectra from the SE and SW cannot be 
well fitted by a single set of n e and T e parameters. The 
emission appears to originate from two distinct dust tem- 
peratures, a warmer component, and a cooler one that is 
more spatially extended. The emission from this cooler 
dust also fills the inside of the SNR (“Inner SNR” re- 
gion in Figure 2). In order to isolate the warmer dust 
component that is only present in the rims, we decided 
to subtract the inner SNR contribution from the SE and 
SW spectra and then fit the spectra with a single tem- 
perature. The final fitting was carried out for regions 
“SE” and “SW’ in Figure 2(b) where region labeled “In- 
ner” was used as the background. This inner region was 
then fitted separately using a background region from 
outside of the SNR shell. The spectra are shown in 
Figures 9 and 10. Prior to fitting, emission lines were 
subtracted and the spectra were extinction corrected us- 
ing the extinction curve from Chiar & Tielens (2006) and 
iVn = 3.5 x 10 22 cm -2 . The best-fit models for the SE 
and SW regions are shown as solid red curves in Fig- 
ure 9, while the best-fit model for the inner SNR region 
is shown in Figure 10. The dotted red curve in Figure 
9 represents the spectrum of the grain size distribution 
prior to any sputtering. The best-fit parameters for these 
models are listed in Table 4, while the x 2 contour levels 
for these parameters are shown in Figure 11, where the 
shaded regions represent values of fl.2.2,4,8,16.32lx min 
X 2 . 

3.5. Gas and Dust Parameters Derived from IR 
Emission 

The diffuse IR emission from the inner SNR region 
seems to originate from a dust population in which the 
larger grains are characterized by a temperature of 55 K, 
while the emission from the SNR shell requires the 55 K 
component in addition to a warm dust component with a 
temperature of 140 K. We used the inner SNR region as a 
background for the SE and SW shell spectra, and in that 
way we have subtracted the contribution from the more 
diffuse, cooler dust component. The residual IR spec- 



Fig. 8. Background-subtracted high-resolution Spitzer IRS LH 
spectra of the southeastern (top) and southwestern (bottom) por- 
tions of the shell in Kes 75. The background was extracted from 
the southernmost LH slit in Figure 3. The observed emission lines 
show no evidence for broadening and most likely originate from the 
background. The emission line intensities are listed in Table 3. 

trum from the warm dust in SE and SW rims requires 
a high plasma temperature and relatively high electron 
density. The best fit value for the electron density for 
both regions is n e = 90 cm -3 , with an electron tempera- 
ture of T e = (1 - -2) x 10 7 K (0.9-1. 7 keV). While the fit 
clearly requires a relatively high density, the temperature 
is somewhat less constrained (see Figure 11). 

The best fit values for the sputtering parameter Aa 
are 0.01 (0.002) pm for the SE (SW) spectrum. These 
values imply that almost all of the smaller PAH grains 
have been destroyed and approximately 50% (25%) of the 
mass of silicate and graphite grains. However, since we 
only have IR spectra longward of 15 pm, this sputtering 
parameter Aa is not well constrained by the model, even 
though the lack of emission at IRAC wavelengths does 
imply that the smaller grains have likely been destroyed 
in these regions. 

While the presence of the cooler plasma component 
is not conclusive from the fits to the X-ray data, it is 
implied by the presence of the 55 K dust component . The 
best fit model to the IR spectrum of the inner SNR region 
implies a lower electron temperature of 1.8 x 10 6 K (0.16 
keV) and a much lower electron density, n e = 2.8 cm -3 . 
The minimum y 2 for the inner region is actually achieved 
for a broader density range of 1-15 cm -3 , but this is still 
several times lower than what is required for the SE and 
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Fig. 9. Dust fitting results to the low-resolution IRS spectra of the SE and SW rims of the SNR shell, with the inner SNR region as 
background. The extraction slits are shown in Figure 2. The spectra have been line subtracted and extinction corrected using the extinction 
curve from Chiar k. Tielens (2006). The solid (dotted) red line shows the best-fit. model with (without) sputtering. The parameters for the 
best-fit models are listed in Table 4. The temperature of the emitting dust in these regions is 140 K. 



Fig. 10. Dust, fitting results to the low-resolution IRS spectra 
of the inner SNR region, shown in Figure 2(b). The spectrum has 
been line subtracted and extinction corrected using the extinction 
curve from Chiar & Tielens (2006). The solid red line shows the 
best-fit model, with parameters listed in Table 4. The temperature 
of the emitting dust in these regions is 55 K. 

SW rims. The Aa parameter is better constrained for 
the inner region. The best-fit value of Aa = 0 implies 
that very little sputtering has occurred in this region (see 
Figure 11). 

4. DISCUSSION 

4.1. Comparison of Derived IR and X-Ray Parameters 

Interestingly, based on the best-fit electron tempera- 
ture to the SE and SW IR spectra, the dust in these 
regions appears to be associated with the higher tem- 
perature (2x10' K) X-ray emitting plasma, rather than 
the 0.2 keV (2xl0 6 K) component as suggested by 
Morton et al. (2007). The association of the IR emis- 
sion with this harder X-ray emission is also implied by 
the spatial correlation of the IR and X-ray emission, as 
seen in Figure 2(b). The IR emission from the inner 
SNR region may be associated with the lower tempera- 


ture plasma component, since the IR and X-ray fits pro- 
vide independent estimates of the n e and T e that are 
consistent with each other. 

The electron density for the two thermal components 
in the X-ray fits was determined from the normaliza- 
tion of the VPSHOCK model, 1 Q~ l4 n e iiuV jiircf 2 cm~ 5 , 
where V is the volume of the emitting region and d is 
the distance to the SNR, set to the value of 10.6 kpc. 
Similar to Morton et al. (2007), we assumed a slab-like 
volume for the X-ray emitting regions. We took the area 
of the SE and SW regions shown in Figure 2(a) and as- 
sumed that the depth along the line of sight is equal to 
the long axis of each region. For an assumed distance 
of 10.6 kpc, the assumed volumes are 15(26) pc 3 for the 
SE(SW) region. We also assume that n e = 1 . 18 jih < as for 
a fully ionized gas. The estimated electron densities for 
the X-ray emitting components are 5.6f~ 1 ^ 2 dfg ^ cm -3 

and 3.lf~ 1/2 dfQ g cm 3 for the ~ 1.5 keV component 
for SE and SW rims, respectively (Table 4), where / 
is the filling factor that accounts for the fact that the 
emitting gas may occupy a smaller fraction of the cho- 
sen volume. The ionization timescale implies a similar 
density, assuming that the shock age is equal to the es- 
timated SNR age of 880 yr. Since the normalization of 
the cooler plasma component is not constrained, we can 
only derive upper limits on the density and gas mass 
in the ~0.2 keV component. The upper limits on the 
density and mass calculated from the normalization of 
the VPSHOCK model (best-fit value plus 1.6tr) for the 
same spatial regions (SE/SW regions in Figure 2(a)) are 

140 dfQ Q 2 cm -3 and 200dfQ g 2 cm -3 , and 43dj’o 2 6 M© and 

110dJo 2 6 M(3> for the SE and SW, respectively. These val- 
ues are listed in Table 4. 

The estimate of n e from the X-ray fits is significantly 
lower than the high density derived from the fits to the 
IR spectrum. If the warm dust component is associated 
with the 1.5 keV plasma, as suggested by implied tem- 
peratures, it may be possible that the X-ray emitting 
gas in the shell is clumpy and concentrated in a smaller 
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TABLE 4 

Dust Model Fits to IRS Spectra 


Parameter n e ( cm 3 ) T e ( K) X 2 Mass (Mq) 


Parameters from IR Fits 

SE shell 

90 

1.8 x 10 7 (1.55 keV) 2.3 

6x 10 -3 (3xl0 -3 ) 

SW shell 

90 

1.0 x 10 7 (0.86 keV) 3.2 

7xl0~ 3 (ox 10 -3 ) 

Inner SNR 

2.8 

1.8 x 10 6 (0.15 keV) 2.6 

0.9 (SE) 
1.1 (SW) 

Parameters from X-Ray Fits (VPSHOCK) 

SE hot 

5.6 

1.7 x 10' (1.46 keV) 

1.7 

SW hot 

3.1 

1.9 x 10 7 (1.60 keV) 

1.7 

SE warm 

<138 

2.3 X 10 6 (0.20 keV) 

<43 

SW warm 

<200 

1.6 x 10 6 (0.14 keV) 

<110 

Parameters from X-Ray Fits (VNEI) 

SE hot 

5.0 

1.9 X 10 7 (1.60 keV) 

1.5 

SW hot 

2.3 

2.2 x 10 7 (2.08 keV) 

1.2 

SE warm 

23 

2.8 X 10 6 (0.24 keV) 

7.0 

SW warm 

<40 

2.3 x 10 6 (0.20 keV) 

<20 

Dust-to-Gas Ratios for the Hot Component 

SE hot 



4xl0“ 3 (2 x 10~ 3 ) 

SW hot 



4xl0" 3 (3 x 10” 3 ) 


Note. — The masses listed for the IR fits are dust masses, while the masses 
listed for the X-ray fits are estimated gas masses. The dust temperatures 
are 140 K (warm component) and 55 K (cool component). The densities 
from the X-ray fits were calculated from the normalization of the VPSHOCK 
model and assuming a filling factor f = 1. Dust masses and dust-to-gas mass 
ratios outside (inside) parentheses represent the best-fit values without (with) 
sputtering. 


volume. If we assume that the density of the clumps is 
90 cm” 3 , as derived from the IR data, we estimate a low 
filling factor that is on the order of a few times 10“ 3 . 
As can be seen from Figure 2(a), the shell does indeed 
appear to be clumpy, so assuming that higher density- 
clumps are present in the shell may not be unreasonable. 
However, the large discrepancy between the X-ray and 
IR derived densities remains a problem. 

4.2. Dust Mass and Dust-to-Gas Mass Ratio 
We estimated the total gas masses in the hot compo- 
nent of the shell to be 1.7 <£(q 2 6 Mq in each of the rims. 
The dust masses were calculated using a grain size dis- 
tribution from the best fit models and the total MIPS 24 
pm fluxes from the SE and SW rims of the shell. The 
fluxes used for the total dust estimates are listed in Ta- 
ble 1. For the warmer dust component we used residual 
fluxes after subtracting the contribution from the inner 
SNR emission. These values are equal to 2.2 Jy (SE) and 
3.4 Jy (SW). We then subtracted these values from the 
total flux in the shell to estimate the contribution of the 
cooler dust component and find these values to be 3.3 Jy 
(SE) and 3.6 Jy (SW). Dust mass estimates are listed in 
Table 4. We find a total dust mass of 6 x 10” 3 df o 6 Mq 
(SE) and 7 x 10” 3 dj 0 6 Mq (SW) for the 140 K dust 

component, and 0.9d% 2 6 MQ (SE) and l.ld% 2 6 M0 (SW) 
for the 55 K component. Since the sputtering amount 
for the hot dust component is not constrained by our 
model, the calculated dust masses assume no sputtering 
and should be treated as upper limits. The best-fit dust 


mass values that include sputtering are listed in Table 4. 
The estimate of the mass of the cooler dust component 
is highly uncertain since we assume that the diffuse IR 
emission all comes from dust with a temperature of 55 K, 
even though this may be true only for the region where 
we have an IR spectrum. A lower (higher) dust temper- 
ature would reduce (increase) the mass associated with 
this component. 

We only estimate the dust-to-gas mass ratio for the 
140 K dust since the evidence suggests that it is associ- 
ated with the hot X-ray plasma. The estimates of the 
dust-to-gas mass ratios are listed in Table 4. Assuming 
that the warm dust component is associated with the hot- 
ter X-ray plasma, we calculate a dust-to-gas mass ratio 
of 4 x lO -3 ^^ 1 ^ 2 f° r an( i SW rims. These values 

are consistent with the dust-to-gas mass ratio of 0.6 1 q) 2% 
that is typically assumed for the Galaxy (Sodroski et al. 
1994). especially considering the uncertainties involved 
in determining the gas and dust masses and the distance 
to Kes 75. If dust sputtering has occurred in the shell, 
the dust-to-gas mass ratio would be a few times lower. 
Since the sputtering amount is not constrained by our 
spectrum, we cannot conclude whether any significant 
dust destruction has occurred in the shell. 

4.3. Origin of the IR and X-Ray Emission 

Morton et al. (2007) attributed the low-temperature 
X-ray thermal component to the forward-shocked 
circumstellar medium (CSM)/ISM, and the high- 
temperature component to the reverse-shocked ejecta. 



